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Introduction: The abundances of noble gases 

found in primitive, carbonaceous meteorites are unex- 
pected when compared with our Sun. Known as Q- 
gases (Q for some unknown carrier dubbed “quintes- 
sence”), this anomaly has remained a mystery since it 
was discovered in 1975. Q-gases are characterized by 
increasing depletions with decreasing atomic number 
(Z) relative to solar noble gases and normalized to 
132 Xe (Figure 1). This Q-gas mass fractionation is un- 
explained, and its investigation is important to under- 
standing the origin of the solar system. 

However, the subject is fraught with controversy - 
in part due to the complex nature of Q and in part due 
to claims of some researchers that cannot be repro- 
duced by other investigators. The topic is discussed in 
numerous places [e.g., 1-4], with models of Q falling 
into two basic categories, both involving carbon en- 
trapment of noble gases. First (Group A), there is the 
conservative two-dimensional view that Q-gases are 
adsorbed or sorbed onto a “labyrinth” of graphite or 
carbon grains [5-9], or they undergo active capture 
onto growing surfaces [6]. Second (Group B), there is 
the view holding to the remarkable property of carbon 
discovered in 1985. Carbon can curl up into closed 
geometries of hexagon- and pentagon-shaped carbon- 
ring configurations, a property ignored completely by 
Group A. Group B thinks of Q as a three-dimensional 
structure of endohedral carbon cages like fullerenes, 
carbon onions, or some class of carbon nanotubes [3, 4, 
10]. Group B does not exclude Group A effects. 

Status: The dismal state of fullerene candidates 

has been reviewed by Heymann el al. [3], while the 
labyrinth group has been equally unsuccessful. Some 
are even saying that fullerenes are irrelevant to cosmo- 
chemistry, a point of view that is certainly false. A 
giant fullerene was discovered in the Allende meteorite 
in 1981 [10, 4, 3] before the advent of fullerene chem- 
istry. If it was not created during or after Allende’ s 
terrestrial impact, fullerenes exist in meteoroids in 
space. 

Following a brief investigation into this circum- 
stance, we report here on a new finding that falls into 
the second group of Q-gas research. 

Q-Gas Abundance: In order to discuss the mass- 
fractionation of Q-gases, their abundance pattern needs 
to be defined. As mentioned, Q contains a noble-gas 
abundance pattern that is anomalous compared with the 
Sun (Figure 1), as well as in comparison with the extra- 
solar (cosmic or interstellar) abundance. Note that Q 


may have two substructures, Q1 and Q2, a suggestion 
dating back to Gros & Anders [11,5]. 

Relevance to Cosmochemistry: Noble gas cos- 
mochemistry must necessarily address the subject of 
endohedral carbon clathrates [12, 13] in view of their 
importance since the discovery of fullerenes and their 
fundamental change in our understanding of carbon 
chemistry. Found originally as the highly stable 
closed-cage configuration C 60 , higher stable fullerene 
allotropes C„ with some having n >100 are known to 
exist in the carbon soot generated by arc-heating of 
graphite. They are created by energetic events, and are 
known to occur naturally and geologically in very low 
concentrations [4]. When produced in the presence of 
an ambient noble gas NG, an endohedral carbon-cage 
clathrate NG@C„ is produced by the reaction NG + C„ 
— ► NG@C„ . Figures representing these clathrates can 
be found elsewhere [4, 12, 13]. 

This theoretical analysis will focus on C 60 per se. 
The carbon-cage clathrates NG@C 60 (NG = He, Ne, 
Ar, Kr, or Xe) are Van der Waals (VdW) molecules 
containing entrapped noble gas NG. Since the noble 
gases are inert and neutral, they interact with matter by 
means of the classical VdW force. The latter, however, 
was not understood until it was derived by London [14] 
as a quantum effect caused by induced multipoles in 
orbital electron density. 



Figure 1. Elemental abundance of Q-gases relative to 
solar abundance and i2 Xe, as a function of mass (from 
[5]) compared to the London or Van der Waals disper- 
sion energy to the fifth power F = E dis 5 for NG@C 60 , 
taken from Table 1. The small dark spots are minimum 
and maximum values for Q from [5]. 


Van der Waals Interactions: When studying the 
behavior of endohedral clathrate systems such as NG + 
C6o NG@C6 o, one tends to begin by thinking in 
terms of VdW radii [e.g., 12] and the Lennard-Jones 
potential U(r) = 4s [(a/r) 12 - ( a/rf ] [13]. The VdW 
contributions (~r 6 ) are attractive at large distances r, 
offset by a repulsive core (~r'~) as r— >0. The - ( a/rf 
term is the London interaction, giving rise to the Lon- 
don force or dispersion force [14]. Because Xe is a 
large atom compared to He, one intuitively would argue 
that Xe@C6o is the most unlikely or unstable clathrate 
in the system - and so on for Kr and Ar. This is oppo- 
site to what is observed in Figure 1 for Q. 

New Developments: Computational quantum chem- 
istry research has made tremendous advances in recent 
years, with much attention paid to endohedral fullerenes 
NG@C6 o- That work includes exchange repulsion and 
electron correlation effects omitted in first-order Har- 
tree-Fock (HF) theory in quantum mechanics. Higher- 
order nonlinear perturbation theory takes the form of 
density functional theory (DFT) and second-order 
Mpller-Plesset theory (MP2) [15], Scuseria [16] and 
Thiel [17] have made some of the first calculations of 
NG@C6 o binding energies derived from DFT and MP2. 

A more recent calculation of London dispersion in- 
teractions for NG@C6 o has been performed by Pyykko 
[18] using MP2. It will be examined here, noting that 
MP2 always lowers the first-order HF state. This re- 
sults in bound states and attractive (negative) London 
terms in Pyykko’ s work. The dispersion energies E dis 
(kJ/mol) are given in Table 1 for Pyykko’s Case 0. 


Table 1 


Noble Gas 

E dis ( kJ/mol) 

Edis 

F = E di f 

He 

-7.5 

0.030303 

2.56E-08 

Ne 

-22.0 

0.151515 

7.99E-05 

Ar 

-71.2 

0.272727 

0.001509 

Kr 

-95.7 

0.636364 

0.104358 

Xe 

-124.9 

1.0 

1.0 


The first thing to observe is that MP2 has lowered 
the dispersion energies E dis as expected, making the 
heavier gases more stable and reversing the trend to be 
compatible with the Q-gas abundance in Figure 1 . Next 
one needs to normalize E dis to Xe, as E dis . Since the 
functional relationship between mass fractionation and 
noble-gas mass in Figure 1 is totally unknown, one can 
ask if there is a straightforward functional dependence 
similar to the Q abundance data. One finds that F = 
E d i S 5 gives such a relationship. This is graphed in Fig- 
ure 1. Except for Ne there is a strong resemblance to 
the Q abundances. 


The same functional relationship used in Table 1 can 
be applied to Thiel’s work on MP2 interaction energies 
(without counterpoise corrections) for NG@C6o bound 
states (-1.0, -5.3, -7.5, -16.1, -17.1 kcal/mol) in [17], 
and one arrives at a similar graph as Figure 1. As 
pointed out by Thiel, the inclusion of a single set of 
polarization functions for the noble gas involved can 
affect the relative energies dramatically. This might be 
relevant to the behavior of He, Ne, and Kr. 

Conclusions: This study does not intend to say that 
the comparison in Figure 1 is a compelling argument 
that fullerenes exist in meteorites or that NG@C f ,n is Q. 
Nevertheless, the resemblance is appealing and de- 
serves further investigation at a much more comprehen- 
sive level, including cases for NG@C„ where n > 60. 
The result is the first of its kind in cosmochemistry, 
derived mathematically from computational quantum 
chemistry. A fifth-order function of the London disper- 
sion energies and clathrate binding energies for 
NG@C6 o bears a strong resemblance to Q-phase abun- 
dances. 
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